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COUPLED-OSCILLATOR  NATURAL  ORBITALS 
Peter  D.  Robinson 

Examples  of  natural  orbitals  in  closed  form  seem  rare  enough  [ 1 ] 
to  justify  this  note,  in  which  such  orbitals  are  exhibited  for  coupled 
oscillators. 

Let  (0  (a,x)}  denote  the  complete  set  of  orthonormal  eigenfunctions 
n 

of  the  oscillator  Hamiltonian 


h(a,x)  = - “ ~ o-^x2,  a > 0,  -oo  < x < oo 

dx 

1 2 

corresponding  to  the  eigenvalues  (n  + — )a  . Specifically 


6n(a,x)  = 


a 


T 3 


it  22  n ! 


H (ax)  exp(-  \ aZx2),  n=  0,1,2. 


where  the  Hermite  polynomials  H^fax)  are  orthogonalized  with  respect 

2 2 

to  a weight  factor  [ 2]  exp(-<*  x ).  Consider  now  the  identical  coupled 

oscillators  described  by  the  Hamiltonian 

„ Id2  1 d2  ^ 1 4.  2 x 2.  1.4  .4..  .2 

H = - 2 TT  - I 7T  + 2 a (X1  + X2)  - 4 (a  - b )(X1  ‘ X2)  ’ 

dxl  dx2 

a > b > 0,  -oo  < < oo,  -oo  < x2  < oo  . 


In  terms  of  the  coordinates 


yl  = ^,Xl  + X2l-  V2  = £(X1  • V ’ 


(1) 


(2) 


(3) 


(4) 
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the  Hamiltonian  (3)  takes  the  separable  form 


H = h(a,yj)  + h(b,y2)  . 


The  ground- state  wave -function  thus  has  a space  part 


1/  \ o / ,Q  ,,  . ab^  . 12  2 12  2 

Wxrx2)  = e0(a,yi)e0(b,y2)  = — exp(-  py,-  jbyz) 


lab  1,2  , 2.  2 2 1.2,2.  , 

= (~|  exp{ - - (a  + b )(x1  + x2>  - -(a  - b jx^}  , 


which  is  the  eigenfunction  of  the  operator  (5)  corresponding  to  the  eigenvalue 

i<*2  + »2>. 

This  function  ijj(xpx2)  is  symmetric  in  x^  and  x.,,  and  is 
normalized  so  that 


S f ^ZdXldx2  = 1 . 


Accordingly,  is  the  Hilbert-Schmidt-type  kernel  of  a compact, 

self-adjoint  integral  operator  on  the  real  Hilbert  space  L^( -<»,<»).  Standard 
integral  equation  theory  [ 3]  tells  us  that  this  operator  has  a complete 
orthonormal  set  of  eigenvectors  {<j>n}  defined  by 

00 

/ ^(x1,x2)<j>n(x2)dx2  = fc-H^Xj),  n = 0,1,2...  , ( 


with  the  eigenvalues  X.  accumulating  at  zero.  Further  the  expansion 

n 

00  00 

-y(xpx2)=/  4'(xp  s)4>(x2,  s)ds  = £ Xn<l>n(x1)4>n(x2) 


u 


is  convergent  for  all  x and  x.  The  functions  4,  (x)  are  precisely 
the  natural  orbitals  for  ^(x^x^),  and  (9)  gives  the  expansion  for  the 
first-order  density  matrix  y(x  x ).  The  "occupation  numbers " \2  must 

i ^ n 

sum  to  unity,  from  (7)  and  (9).  The  expansion  for  ^(x^x^)  itself,  namely 


00 


must  also  be  convergent,  at  least  in  the  L 2 norm  sense. 

In  the  case  of  the  function  given  in  (6),  the  expansion  (10)  can  be 
found  directly,  and  is  convergent  for  all  Xj  and  x^  Mehler's  formula  [ 2] 
states  that 


exp  ( - 


r7(t‘ + ‘^+ r?  *>l*J = (1  - z2,i  n?„  fz H- wv> 

-1  < Z < 1.  -00  < t < 00  -00  < t <00 

’ 1 ’ 2 


00 


(ID 


Setting  t^  = ax^,  t^  = ax  and  multiplying  through  by  exp{-  - <*2(x2  + x2)}, 

2 1 2 ’ 

we  obtain  from  (2): 


a f II  + z 
— exp(  - 


'Hi-,2 


2l  «2  2 2 <j2z 

T<xi  4 x2)  + 2r"ixix2 


Z ~ n 

(1-z  )2  Yj  zn0  (ar,x  )0  (ar,x  ).  (12) 

a •*  An  l 


1 - z 


n=0 


With  the  choices 


a = (ab)2,  z 


a - b 
a + b 


(13) 


the  left-hand  side  of  (12)  becomes  4j(x^,x2)  as  given  in  (6).  Thus 


-3- 


* 


i 

' 

I 


I 

I*. 


r : 


*<*,.«*>  ■ l ftr?  <-«"( - 

n=0  ^ 


so  that  for  the  ground- state  of  the  coupled  oscillators  we  have  the 
natural  orbitals 


<b  (x)  = 0 ( \/ab,x)  = 
n n 


\rab 


Ltr  z2  n U 


Hn(NJabx)exp(-  abx2)  , 


and  corresponding  eigenvalues 


2\/ab  . n/ a - b \" 

' ' la  + b)  ‘ 


n a + b 


\ = 
n 

2 

The  orbitals  (15),  and  the  occupation  numbers  Xn,  can  also  be  obtained 


1 


by  working  from  the  expression  for  - y(x, ,x  ),  which  is 


2 

4 .4  , 2,  2, 


2 . 2,2 


1 . . fir  . 2 . . 2.1  2 J (a  + b + 6a  b ) , 2 . 2,  (a  - b ) 

iV<x1,x2)=ab[-(a  tb  )J  exp(-  (X1  + *2»  * „ 2 . 2 *1 


8(a  + b ) 4(a  + b ) 

Extension  to  three-dimensional  oscillators,  for  which  Davidson  [ l] 
has  discussed  the  first-order  density  matrix,  is  straightforward. 
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